On the Equivalence of Quantum and Classical Coherence in Electronic Energy 

Transfer 
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To investigate the effect of quantum coherence on electronic energy transfer, which is the subject 
of current interest in photosynthesis, we solve the problem of transport for the simplest model 
of an aggregate of monomers interacting through dipole-dipole forces using both quantum and 
classical dynamics. We conclude that for realistic coupling strengths quantum and classical coherent 
transport are identical. This is demonstrated by numerical calculations for a linear chain and for 
the photosynthetic Fenna-Matthews-Olson (FMO) complex 
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INTRODUCTION 



The role of electronic excitation transfer (EET) in the 
photosynthetic process has been emphasised since the 
1938 study by Franck and Teller Already in that 
paper EET was suggested to occur through the coher- 
ent quantum motion of Frenkel exciton waves [ij. Re- 
cently renewed interest j3l-[lo| has been awakened by 
spectroscopic observations pointing to long-lived quan- 
tum coherences in the process of energy transfer in pho- 
tosynthesis [llj. In particular it has been claimed that 
the speed of EET is enhanced by quantum effects [loj . 
which themselves are manifestations of quantum entan- 
glement. The analogy has been made to advantages in 
quantum over classical information processing. Consid- 
erable publicity has been accorded these claims which 
have been suggested to be the first example of 'quantum 
weirdness in physiology' 12| . Hence it appears timely to 
test this hypothesis on the simplest exciton model of an 
aggregate of dipole-coupled monomers, which we treat 
by both quantum and classical mechanics. We demon- 
strate, both numerically and analytically, that for dipole 
coupling strengths relevant to EET in molecular aggre- 
gates, the coherences in quantum transport (from the 
Schrodinger equation) are identical to those occurring in 
classical transport according to Newton's equation. Cor- 
respondingly, in this fundamental example, the effect on 
the EET mechanism arising due to quantum entangle- 
ment is mirrored by the coherence of classical dipole- 
dipole EET. 



Although coupling to internal nuclear degrees of free- 
dom and to the environment plays an important role in 
EET (see e.g. Refs. [i,li,[i3-[ii), here we will concentrate 
on the most basic situation in which quantum coherence 
appears. Accordingly we consider the idealised case of an 
aggregate of monomers, each having only one electronic 
excited state, whose purely electronic degrees of freedom 
are coupled by dipole forces, shorn of any complications 
arising from dissipation or decoherence due to vibronic or 
environment coupling. Quantum transport will be deter- 



mined by solving the time-dependent Schrodinger equa- 
tion (TDSE). 

To treat the quantum system classically we must 
relate the classical motion to quantum properties. Here 
we consider the monomers to consist of an oscillating 
classical electron whose dipole strength will be set 
proportional to the quantum oscillator strength. Exactly 
the model we use here has also been used to discuss 
excitonic light absorption by solids and liquids (see e.g. 
Refs. 



20l-l22|) and e nerg y transfer on dimers 



nanoparticle arrays [25 



and 



QUANTUM MECHANICS 

The time-dependent Schrodinger equation (TDSE) for 
a Frenkel exciton on an aggregate is 

zft||*(i)) = (Ho + V)|vl/(t)) (1) 

where Hg is the sum of the hamiltonians of non- 
interacting monomers and V is the total potential energy 
of the pairwise interactions between monomers. Since we 
consider the propagation of a single electronic excitation 
along the aggregate, we expand the full time-dependent 
wavefunction as, 



(2) 



where e„ is the single-monomer transition energy and 
1 7r„ ) denotes a state in which monomer n is electronically 
excited and all other monomers are in their ground state. 
The full aggregate ground-state energy is set to zero. The 
matrix elements of Hq are taken to be 



( 7r„ |Ho| 7r„ 



(3) 



To connect with the classical model of oscillating dipolcs, 
we will take the interaction operator V to be the dipole- 
dipole interaction 



^nrn — ( |^| '^m ) l^n^n -^nm -^rnl^m 



(4) 
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where /i„ and e„ denote the magnitude and orientation 
respectively of the transition dipolc matrix element of 
monomer n and 



7nn ' 



(5) 



where Rmn is the direction of the vector separation, with 
length Rmn, of monomers m and n. Substitution of 
Eq. (12]) in the TDSE leads to a set of coupled equations 
of first order in the time derivative, 



(6) 



These are the quantum equations. 

CLASSICAL MECHANICS 

In the classical case we consider that the coupled quan- 
tum transition dipoles are modelled by classical oscilla- 
tors in the same geometry. The Hamilton equations of 
motion for linearly- interacting oscillators of mass and 
frequency a;„ are, 

Pri /_N 
p„ = - Mnl^lxn - ^ /•i'nmaJm (8) 



where the Knm sue coupling coefficients determined by 

Xn{t)K^7nX777.{t) = flfit) ■ T77777 ' A^n^ (9) 

and the r.h.s. has to be related to the quantum mechan- 
ical dipole-dipole interaction. 

The treatment of quantum transition dipoles, cor- 
responding to electronic excitation and de-excitation, 
as classical electric dipoles requires incorporating the 
monomer quantum properties into a classical model. 
Here we follow Fano |2(| in defining classical dipoles 
whose dipole moment corresponds to the quantum tran- 
sition dipole moment. To this end we first introduce the 
dimensionless quantum oscillator strength 



(10) 



where rrie is the oscillator (electron) mass, e is the elec- 
tron charge and e„ = hiUn relates the transition energy 
to the monomer oscillation frequency. Eq. (fTO]) allows us 
to represent the classical dipole moment in terms of the 
quantum transition dipole moment fin- That is we take 
the classical dipole moment of monomer n to be given by 



class 



(11) 



where e„ is the direction of the transition dipole of the n- 
th monomer. After identifying M„ with rrie the classical 
coupling strengths are then given by. 



■T, 



(e^V/n/m)™e 



(12) 



These coupling coefficients can now be related to the 
quantum mechanical Vnm, given in Eq. (|4]). By inserting 
Eq. (dni) into Eq. dH]) on finds 



2V„. 



(13) 



To connect to the quantum equations, we introduce the 
dimensionless complex amplitude 



Zn[t) = Xn{t) + ipnit) 



with 



2h 
1 



Pn 



(14) 

(15) 
(16) 



One then obtains 



iZn = UJnZn + ^ ^^'" Re(Zr») (17) 



This equation has to be compared to the quantum equa- 
tion ([5]). The difference is a factor of two and the ap- 
pearance of Ile{zm) instead of Zm in the coupling terms. 
The resulting difference in the dynamical equations is 
best seen by first taking the time derivative of Eq. (fT7|) 
to obtain, using p„ — Xn/iOn, 



Z71 — —Ldf^Zn — UJ,i 



(18) 

By contrast, taking the time derivative of the TDSE 
Eq. we obtain a second order equation for the quan- 
tum amplitudes. 



2 2V^77i 



— OJn) — T~^7i 



E 



(19) 



Scaling the time as a;„f this equation is, 



E 



c„ 



(20) 



^71 ^71 



We note that the last term in this equation is of sec- 
ond order in the normalized coupling strength Vnm/{£7%) 
and can be neglected when this coupling is small. This 
we call the realistic coupling approximation (RCA) . The 
penultimate term is zero for identical monomers em = 
and small when the fluctuations in the transition energies 
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are small relative to the mean. Under the same condi- 
tions the last term in the classical Eq. ([T8)) is small. In 
these two approximations, we can drop the last two cou- 
pling terms in Eq. and the last term in Eq. ([T5)) . so 
that these two equations become identical. To this order 
then the classical Eq. ([T7| and the quantum Eq. © are 
equivalent and we can associate, up to a normalization 
factor, the quantum amplitudes c„ with their classical 
counterparts z„. The classical "probability of occupa- 
tion" P,f '''''' (t) is then defined as 



(a) 



(b) 



class 



it) 



(21) 



In realistic cases of aggregates involving monomer 
molecules one has the situation that the interaction 
splitting is much less than the transition energy i.e. 
Vnm <C e„ V n,m, i.e. the RCA is valid. Indeed, 
in the light harvesting complexes the electronic tran- 
sition energy of the chlorophyll monomers is ^ 12 000 
cm~^, whilst the dipole-dipole interaction between the 
chlorophylls is in the order of a few hundred cm~^ 
(27I I. Hence we can conclude that all coherences of 
propagation in the quantum case are also contained in 
the classical propagation for realistic coupling strengths. 
Note that the quantum equation Eq. ^ can always 
be simulated by classical oscillators coupled linearly in 
position (or dipole strength), so long as all {Vnm/^n) 
are small and the relative fluctuations in the e„ are 
small. We emphasize that this is quite different from 
a standard mapping [H, [2^ of electronic degrees of 
freedom to classical oscillators. This mapping gives 
an unrealistic complex classical Hamiltonian involving 
additional momentum couplings. Only in RCA do 
the equations (fTT]) for coupled classical electric dipoles 
reproduce the quantum amplitudes of coupled transition 
dipoles i.e. only in this approximation is Eq. (jl7p a good 
approximation to Eq. ([6|). This is the major result of 
this analysis. 



COMPARISON OF QUANTUM AND CLASSICAL 
MECHANICS 



In the following we will demonstrate the applicability 
of the realistic coupling equation with two examples. The 
first one is a linear chain of identical monomers and iden- 
tical coupling between neighboring monomers. Writing 
Kim = V and taking only nearest neighbor coupling into 
account the equation ([6|) reduces to 



.dcn 



1 



(c„_i(t) + c„+i(r)). 



(22) 



Here the dimensionless time r = {2V/K)t has been in- 
troduced. One can show that the solution of these equa- 
tions with initial condition c„(0) = 5no are just the Bcssel 
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FIG. 1: The probability as a function of time in units of h/2V 
of excitation of monomer n for a linear chain. Blue: quantum 
result, red: classical result. Left figures for V/e = 1/40 and 
right figures for V/e =1/6. 



functions c„ = J„(2V^i/fi,) exp(m7r/2). Then, beginning 
on monomer zero, the excitation probability of monomer 
n is 



Pn{t) = Jl{2Vt/h) 

This is the result obtained e.g. in Refs. [13, Hc 
the analytic result 



^Y.^''jl{2Vt/h) = {2V/hft^ 



(23) 



From 



(24) 



one infers a constant "velocity" of electronic propaga- 
tion given by {n?) /t = 2V/h. Hence in this sim- 
ple model the transport of electronic energy is fully co- 
herent. Note that, in nearest-neighbour coupling, the 
scaled time h/ {2V) is the natural time unit since it is 
just the mean time of EET from one monomer to its 
neighbours. A detailed study of this coherent quantum 
transfer is given in Ref. [17[. Within the RCA the clas- 



sical equations predict identical coherence and constant 
velocity of propagation. We have tested the validity of 
the RCA by performing exact numerical evaluation of 
the full classical equations Eq. ([7]) and Eq. ([8|) for this 
case of identical monomers and nearest-neighbour inter- 
action. The P^^'^^^ (r) for the classical energy propagation 
along a chain of 19 monomers are compared to the quan- 
tum Bessel function prediction in Fig. [T] where results 
for two different coupling strengths are shown. One sees 
(Fig. [I}i) that for a small value 1/40 of the parameter 
V/e, nevertheless one that is typical for photosynthesis. 
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the exact classical probability that energy resides on a 
given monomer at a given time is indistinguishable from 
the quantum Bessel function result. Even for Vje = 1/6, 
which is unrealistically large, one sees (Fig. [T|d) that the 
classical result still follows closely the envelope of the 
quantum result, although some deviations are beginning 
to appear. Hence, for realistic coupling, these numerical 
results fully justify our analytic approximation leading to 
identity of classical and quantum EET. 

In the previous example of a linear chain all the 
nearest-neighbour coupling strengths and transition en- 
ergies have been taken to be equal. To demonstrate that 
the classical and quantum results also coincide in more 
general cases, we have calculated energy transfer dynam- 
ics for an aggregate where the coupling elements and 
transition energies cor resp ond to those usually assumed 
for the FMO complex [31[ . Coupling to the environment 
is neglected. The complex consists of seven monomers 
and we have calculated the probability that excitation, 
starting localised on one monomer, has reached other 
monomers at later times. Typical results are shown in 
Fig. [5] for excitation initially localized on monomer 1. 
Since monomer 1 couples strongly to monomer 2 only 
the probabilities for these two monomers are of the same 
order. Also shown is the much lower probability of ex- 
citing monomer 3. The probabilities of excitation of 
other monomers are not shown but the agreement be- 
tween quantum and RCA results is equally good. The 
left column shows results using realistic FMO couplings 
and transition energies quoted in Ref. [sH. One sees an 
oscillatory motion of excitation but most important, the 
classical and quantum results are indistinguishable and 
appear as a single curve. In the right column the tran- 
sition energies have been artificially reduced by 12000 
cm~^ so that they become similar to the coupling ener- 
gies. Then the RCA is no more valid and classical results 
(red) arc not identical to quantum (blue), although even 
in this case the overall agreement is still very good. Simi- 
lar agreement has been found for all other monomers and 
also for different initial conditions. 

What is important here is not to infer that, since the 
quantum result agrees with the classical, there is no quan- 
tum entanglement. Rather it is the classical amplitudes 
and phases which mimic this entanglement. For the one- 
exciton manifold of pure states considered here, entan- 
glement is expressed simply by the possibility to write 
the wavefunction, or parts of it, in a single product form 
or not. A quantitative measure of the entanglement 
between sites i and j is the concurrence [s^ which in 
the present situation is simply proportional to the ma- 
trix elements of the corresponding density matrix, i.e. 
I/Ojjl = |c*Cj|. As example, for the linear chain we have 
the analytic result \p^j{t)\ = \Ji{'lVtlh)Jj(2Vtlh)\ de- 
scribing the development of entanglement in time. Note, 
however that the classical complex amplitudes will lead 
in RCA to exactly the same result. Indeed, quite gen- 
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FIG. 2: Time dependent electronic occupation of monomer 1 
(top), monomer 2 (middle) and monomer 3 (bottom) in the 
case of the FMO parameters. Initially the excitation is lo- 
calised on monomer 1. The right column shows results where 
the transition energy is arbitrarily decreased by 12000 cm~^. 

erally wc find that the agreement between quantum and 
RCA classical complex amplitudes, which determines the 
coherence, is equally close as for the population compar- 
isons shown in the figures. 

To summarise, the question as to whether quantum 
mechanics has an effect on EET on molecular aggregates 
appears from the results presented here to be answered 
with a definite 'no'. One can conclude that the entangle- 
ment in the quantum case leads to no essential change 
from the coherence characteristics of classical electronic 
energy transfer. The results presented above are ob- 
tained by neglect of the influence of an environment. 
We have also performed calculations including dephas- 
ing processes in the quantum and classical equations and 
found similar good agreement. Also, in a recent paper, 
Zimanyi and Silbey [33[ have studied the related prob- 
lem of energy transport between just two monomers of 
differing character in the presence of environmental ef- 
fects. They show also in this case that a purely classical 
model leads to the same energy transfer probability as 
the quantum result. We conclude that some care must 
be taken when attributing coherence in energy transfer 
as due to quantum mechanical effects. 
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